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We study the Regge trajectories and the quark-antiquark energy in excited hadrons 
composed by different dynamical mass constituents via the gauge/string correspondence. 
First we exemplify the procedure in a supersymmetric system, D3-D7, in the extremal 
case. Afterwards we discuss the model dual to large- A'^c QCD, D4-D6 system. In the latter 
case we find the field theory expected gross features of vector like theories: the spectrum 
resembles that of heavy quarkonia and the Chew-Frautschi plot of the singlet and first 
excited states is in qualitative agreement with those of lattice QCD. We stress the salient 
points of including different constituents masses. 
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1. Motivation 

One of the main features of vector-like tlieories is tlie self-interaction of the spin one 
gauge fields. With the inclusion of fundamental matter, colour charged, two competing 
contributions appears to the quark antiquark potential: one, attractive, is due to the 
interchange of the gauge bosons while the second, repulsive, comes from the self-interaction 
of these gauge bosons. The net effect is the formation of an electric flux between the two 
colour sources. The confined, electric flux tube connecting massless quarks can be modelled 
by a string, which represents one of the simplest hadron models When the string 

tension is constant it leads to a simple relation between the energy and the spin of the 
string, E ~ i.e. Regge. From the quantum states assignment, Regge trajectories 
are based on flavour independence arguments and exchange degeneracy, i.e. J^'^ and 
(J + states belong to a common Regge trajectory. Even though it is well known 

that more realistic Regge behaviours include a positive curvature and a nonzero intercept 

i- 

It is our aim to recast the features of these models in the gravity /field theory duality 
framework. For that purpose we shall deal with the meson Regge trajectory in the light 
of holographic duals of confining supergravity backgrounds. In the first holographic setup 
1^,0 the addition of matter was achieved by higgsing the initial gauge group: A stack 
of D3-branes was considered at the origin giving rise to a gauge group SU{N). Then 
one of the branes was pulled out up to the boundary breaking the initial gauge group to 
SU{N — 1) X U{1). The string length stretching from the origin up to the boundary is 
proportional to the mass of the quark attached at one of the string ends. Then the final 
picture was an infinite heavy quark at the boundary. Considering two of such strings a 
more stable configuration, merging the lower end points, is possible giving a bound state 
of two infinite masses. 

The addition of finite fiavoured matter is achieved by means of an open string sector 
0. The end points of these are supported by Nf nonbackreacting filling space-time D- 
branes, Nf <C Nc —>■ oo . In most of the cases, D-branes carry RR-gauge fields and hence a 
nonzero charge. A possible mechanism to stabilise them at finite distance is by wrapping 
trivial cycles [^H. Then, without creating instabilities, negative massive modes prevent 
the slipping of the D-brane out of the wrapped trivial cycle. The final result are a set of 
mesons (in the large- A'c limit) consisting of bound state poles. Neglecting meson widths 
correspond to disregard mixing between qq and other sectors of the theory as glueballs or 
qqqq. 

The meson spectrum can be divided in two different sectors: i) A sector with small 
(J = 0,1,2) angular momentum, that can be tackled considering small fluctuations of 
the probe-brane. ii) Stringy states at large-J corresponding to excited meson states. In 
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the latter case they are represented by rotating Wilson-loops. Notice that this inherently 
implies a semi-classical limit. It is then natural to assume that isolated hadrons exist at 
sufficient energy to be treated semi-classically p. 

We shall consider the spectrum of the stringy states in a generic case: an hadron 
formed by two different constituents masses. This will allow us to study some features of 
the hadron excited spectra as a function of the masses in a more realistic set up. There is, 
for instance, little knowledge of heavy-light charmed mesons. While at low-energy there is 
no significant change with respect to the degenerate case, at high-energy the entire profile 
changes. In the cases where the metric has been obtained via compactification we argue 
that in the vicinity of the high-energy region there is the opening of the compactified 
dimensions, and hence the dual field theory is changed. 

In order to introduce notation and the main work line we deal with an extremal 
supersymmetric model, D3-D7. We discuss the spectrum, the energy-momenta relation and 
the quark-antiquark energy. We pay especial attention to the role of mass subtraction. We 
turn then to evaluate the above features in a nonsupersymmetric model, D4-D6, comparing 
their qualitative behaviour with the experimental ones. In addition we comment on the 
quark-antiquark potential screening and on the 't Hooft line. 

Even if we treated M = 2 and A/" = theories it would also be interesting to 
perform a similar analysis for M = 1 theories although, being confining, it is natural to 
expect a similar behaviour to that of the non-supersymmetric case. The embedding of the 



corresponding flavour branes was studied |lC]in the case of a non-conformally flat setup 



and in [0,|T2|,|T^ for the theories on the conifold. 



2. The D3 - Dl brane model 

We start with a supersymmetric system that allows a fairly good analytic control over 
the full calculation. To obtain a system with matter in the bi-fundamental representation 
we make use of the observation that a holographic dual to 4-dimensional super- Yang-Mills 
(SYM) theory can be obtained in the near-horizon limit of the D3-D7 system 

The set up consists in a stack of Nc D3-branes at the origin expanding in four- 
directions (0, 1, 2, 3) and warping along the transverse ones. This give rise to an A/" = 4 
SYM fleld theory. Another set of Nf D7-branes are located along (0, 1, 2, 3, 4, 5, 6, 7) and 
warping in the rest of dimensions. The insertion of this last stack of branes breaks the ini- 
tial A/" = 4 symmetry to A/" = 2 SYM at low-energy. We shall consider the case Nf <^ 
then the D7-brane backreaction can be ignored. 

The D3-brane system can be obtained from type-IIB string theory on the AdS^ x 
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background 

ds' = fiU)-'/^-dxl + dx')+f{Uy/'dy\ f{U) = '^, (2.1) 

with X = {xi,X2,X3) ,y = {x^, . . .xg) , p'^ = and TZ^ = Qs Ned . The induced 
metric on the D7-brane embedding is that of AdS^ x 5""^ , i.e. the D7 fiUs aU the AdS space 
but is only a contractible cycle on S^. This breaks the initial SO [6) symmetry on the 
5"^ to SO {4) X SO{2). The former corresponds to the isometry group of the S^ while the 
latter to the rotations on the (xg, xq) plane. 

The fundamental massive gauge theory is found by separating, a distance £, the D7- 
branes in the xs — xq plane, x| + x| = £^ 

ds^ = ^''^f [-dxl + di") + ^y^^ dOi , dQ\ = {dXJ'^ + t/2d02) . (2.2) 

In fiat Minkowsky space the quark mass is proportional to I. At the ultraviolet region, 
V ( p.2|) asymptotes to AdS^ x S^ providing a conformal invariant gauge theory. In 
reducing the energy, decreasing t/, conformal invariance is lost. 



2.1. Spinning mesons 

We start by determining the possible spectrum for open rotating strings with ends 
attached to the D7-branes t). The string worldsheet action is parametrised by (r, a) as 

S = [ drdaC (X, d^X, d^X) . (2.3) 



The string configuration we choose is the same as the one described in |T^] but with the 
insertion of a second brane probe. It will lie at a fixed point in 5"^, rotating with constant 
angular velocity uj within a 2-plane in E'^ and extending in the direction z = TZ'^ /U. The 
relevant part of the metric is then 

ds^ = % {-de + rfr2 + r'^de^ + dz^) . (2.4) 
z^ 

Defining Zm^ as the position of the flavour brane in terms of the new coordinate z, the 
mass of the quark is 

m, = — — . (2.5) 



2'nzr 



The chiral symmetry breaking of the model was studied in [14|. 
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The rotating string describes a cigar-like surface with the configuration 



(2.6) 



Under the latter assumption, ( p.6|) , the Lagrangian density takes the form 

7^2 1 



(2.7) 



where we have already substituted u) for 6 and denoted by primes. The equations of 
motion for z{a) and r(a) are equal and they reduce to 
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^'2 _|_ j,'2 



f I II II i\ , i 

72 (^r - z r ) + UJ 



2- = 0. 



(2i 



In addition to ( |27q ) one has to supplement the system with the proper set of boundary 
conditions (b.c). In order to find the b.c. to be applied we use the minimum action 
principle: the string must be stationary under small deformations of the world-line 



SS = -2 J dT da Sz ^ V(k2 - u;2^2) (^/2 + ^/2) _ ^2 y 



'2 — ci;2r2 



, , ^ r' I k"^ — uP't'^ 
-I- / drdr — 



z" 



r'2 + ^'2 



r/2 



, , J. z' / ^2 — c<;2r2 
+ drdz ^ 



^/2 I ^"^2 



r/2 



0. 



(2.9) 

/ Z" V -i- Z'~ 

cr = -7r/2 

Each end of the string lies on one of the D7 branes, i.e. ^-coordinates obey Dirichlet b.c. 
While for 



r obeys a Neumann b.c. . From (|2.9| ) and in the case of subluminal 
velocities this implies that the string ends hit the D7-brane orthogonally. 
The e.o.m. are easily worked out by fixing the gauge 



i) r{a) = a 



ii) z{a) 



12 



J^a , 2 

+ - = . 



l-a;2^2 



r 

2- + 



z 



r = 0. 



^2.10) 



:2.111 



1 + r'2 (J 1 — (^2j,2 

Notice that (|2.10|) , ( |2.11j ) are identical to the ones given in [|15|. Our only concert here 
has to do with the b.c. imposed on these 



Zmq , = fixed , daZ\a=±^ -* OO 



(2.12) 



In the limiting case, Zm^ = zm^ we recover the results presented in []T5 
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The system (|2.7P does not depend explicitly neither on the time coordinate nor on 6. 
This implies the existence of two conserved quantities. Moreover using reparametrisation 
invariance, k = 1 , the energy and the spin are given by 



■^=7^377/^^(7) 1/^-272 • (2-14) 



□ □ □ 



.25 0.5 0.75 1 1.25 1.5 1.75 



Fig. 1: Chew-Frautschi plot for niq = 1 and Mq = 1,2,3. The dashed lines 
represent the limits (2.18) and (2.29) for the equal and degenerate case. We also 
have plotted (horizontal thick line) the energy of a free system of quarks of energy 



2mq . 



We shall study the behaviour of ( pTO|) , ( pTT|) , ( pT3D and (|]l|). The most straight 



forward way is a numerical analysis. In fig. 1 we show the energy versus the squared root 
of the angular momentum b. As one can see the degenerate system behaves at low-angular 
momenta in a completely different manner as the different masses case. Even though the 
analysis is completed it would be illuminating if we recover analytically the result in two 
limiting cases. We follow closely [|T3[ by slightly adapting it to the present case. 



2.1.1. J < VgsNc Mesons 

Let us analyse the J limit. As depicted in fig. 2, the configuration corresponding 
to very small J is a string that directly stretches in the z direction between both fiavour 
branes, without bending. Obviously, this is qualitatively different of what happens in the 
single flavour case. In order to have J — 0, we need r ^ (for all the string worldsheet), 



b Even if we plot dots, this by no means is because any quantisation of the angular momenta. 
Those intent only to be representative points. 
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Fig. 2: The profile of the string for J <C \^gsNc for ruq = 1, Mq = 2 and 
a; = 4. The dashed lines represent the analytic continuation of the string profile. 
As suggested in ||l5|| those lines correspond to gluons shells. Contrary to the case 
presented in the above reference in the case of nondegenerate mass these gluons 
shells are not longer symmetric. The horizontal lines represent the probe-branes. 



as one can realize after a glance at ( p.l4|) . Therefore, we are in the limit corresponding to 
rotating short strings. Neglecting higher order terms in r, (|2.11|) becomes 

r" -2- +uj^r = 0, (2.15) 
z 

with r{zmq) = e and where only the leading order in epsilon is taken into account. 

Supplementing this equation with the boundary conditions r'^Zm^) = 0, r{zmq) = e 
we find 

r(z) = [zuj cos {{zm^ - z)uj) + sin {{zm^ - z)uj)] . (2.16) 

Moreover, by fixing r'{zMq) = 

uj = . (2.17) 

This is a limiting upper value of the angular velocity, contrariwise to the one-mass case, 
where uj ^ oo when J ^ 0. 

Notice that niq r{zmq) = Mq r{zMq) ■ This is the momentum conservation relation for 
two rotating masses, i.e. the string does not carry any energy in this limit acting only as 
a geometrical constraint. 

The energy-spin relation is given in this case by 

^ , , 27r^Q;' rriaMr, 

E^Mq-mq+ ^ J. 2.18 

IZ^ Mq -f rUq 

This linear behaviour relating the energy and the angular momentum at high energy has 
no field theory interpretation. 



2.1.2. J > ^JqsNc Mesons 
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Fig. 3: The z{r) profile for the degenerate, l.h.s., and nonequal mass case, r.h.s. . 
The dashed lines represent the analytic continuation of the profile following eqs. 
( 2.101 ), ( 2.11 ). The horizontal solid lines represent the two flavour branes. z = 
represents the AdS boundary. 



The other possible limit where ( p.l3| ) and ( |2.14| ) can be tackled analytically is w — 
which corresponds to rotating long strings that give large J mesons. We shall generalize 
the analysis of in order to incorporate the effect of having two different quark masses. 
For small angular velocity, the profile of the string must be similar to that of a static string 
hanging between both branes. This is given by 



rst{z) =ro + 



20 



dx 



X 



(2.19) 



where 2:0 ~^ oo in the limit we are considering and it refers to the extremal point. We have 
defined ro as the r coordinate of the point where the string bends back. It cannot be set 
to zero as r = is by definition the center of rotation of the string and both points do not 
coincide for nondegenerate masses (see the asymmetric plot in fig. 3). Let us point out 
that the distance between the string edges is 



with 



C 



'° dy v^r (3/4) 



(2.20) 



(2.21) 



'1 y^y^rrrr r(i/4) 

By considering r{z) = rst{z)+Sr{z) and keeping the leading order in Sr{z), uP' , ( |2.11| ) reads 



z 1 — 00 r^^^ 



(2.22) 



This equation can be straightforwardly integrated 

-v2 



Sr'{z) 



{zl-z^)i 



const -\- uj'^Zq 



20 



x^ 



1 — uj'^r'ii 



-dx 



St 



X) 



(2.23) 
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and in addition has to be supplemented with the boundary conditions (|2.12|) . They read 

r'^,{z^^) + 5r' {zm,) = K.izM^ + Sr' (zmJ =0. (2.24) 

As we are in the long string regime, one can approximately think on the profile of the 
string in the following way: It stretches vertically from ZMq to Zq at constant tm^, then 
stretches horizontally from to at constant zq and finally stretches again vertically 
from Zq to Zruq at constant . A glance at ( |2.23| ) reveals that the integral is dominated 
by the vertical regions i.e. x <^ zq. Moreover, the constant can be set to zero as the point 
with 6r' = is in the horizontal region, at least when the quark masses are of the same 
order of magnitude (this reasoning fails when ^ oo, because when there is an infinite 
mass r!.ti^Mg) = so dr'{zMq) = 0, i.e. the Sr' = point is in the vertical region). 
Then, one finds from ( |2.24| ) 

^2 ^ ^ ^ JM^ ^ (2.25) 

where we have used the relation zq ^ Zm^, ZMq and, as can be proved from ( |2.2C1| ),( |2^^ ), 
iv'^rlj^jiv'^r'^^ <^ 1. From (|2.25| ) it is immediate that 

mqrmq = MqTMg ■ (2.26) 

This is the momentum conservation relation for two rotating nonrelativistic masses, and 
signals that the horizontal part of the string carries a neglectable amount of energy. 

Using the relations above, it is easy to compute the angular momentum from (|2.14|) to 
leading order 

J^J1^ + .... (2.27) 

na Zq uj 



The energy (|2.13|) , considering the leading and subleading orders is 



7^2 / 1 1 C , , , 

E^- — , + +.... (2.28) 



27ta' \Zm.g ZMg Zo 

Using (|2.20|) , (|2.25| ) and ( |2.27| ) this expression can be recast as 



E = + mg - ^binding , ^binding = + M ) 4J2 ' "^^'^^^ 



where 



We deal with a system of two nonrelativistic particles bounded in a Coulombic potential, 
V{r) ~ — K^/r , where only a finite energy cost is needed to separate the two masses. Notice 
that the binding energy is proportional to the reduced mass of the system. At relatively 
large separations, J ^ , the hadron mass asymptotes to Mg + ruq . While for given masses 
and strong coupling -^binding ^ {Mq + niq), hence been the string stable. When the limit 
Mq —>■ oo with niq held fixed is taken the binding energy must become independent of the 
coupling Ijl^ and our procedure fails. 
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2.2. Quark- antiquark energy 

The Wilson loop provides an order parameter for pure gluo dynamics. Even if in the 
presence of dynamical quarks this feature is washed out, it still shows a strong sensitivity 
in the transition region and therefore is widely used to describe crossover in full QCD. A 
possible gravity analog for infinitely heavy-quarks was presented in 

The addition of dynamical flavour test sources , D7-branes, allows for the inclusion 
of light-quark masses . It is of interest to investigate the behavior of the qq energy as 
a function of their relative separation in this case. One of the expected new features is 
the screening in the potential due to light-quarks. To compute the qq energy we locate an 
open string with the two ends on the probe-branes and stretching from them to the origin 
in the U coordinate. At some point Uq, before hitting the horizon, it bends back. The 
relevant part of the metric reduces then to 



r^2 ^2 

ds'^ = a' — (~dt'^ + dr^) + -r-rdU' 
IZ U 

The string configuration is parametrised as 

U = U{a), r = a, t = T. 
This leads to the lagrangian density 



(2.31) 



1 rrms.^ / 7-7-4 

C = / dr\ t/'2 + ^ . (2.32) 



The independence of (|2.32 ) w.r.t. the world-sheet variable a implies the existence of a 
conserved quantity 

= n^U^ . (2.33) 



u 




Fig. 4: The profile of the string in the U — r plane for different lengths. The figure 
on the left corresponds to long strings and the one on the right to short strings. In 
both cases we take niq = 1, Mq = 2 
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In fig. 4 we have plotted the profile of the string obtained by numerically solving ( |2.33 ). 
For long strings there is a turning point while short strings stretch directly between probe 
branes. Notice that the long string limit is achieved by approaching Uq . 

The separation between the two ends of the string can be easily computed from ( |2.33| ) 

L = 2Tr/ o r^ + TT / ^ o ■ (2-34) 



Uo J I y^^y^ - 1 Uo J2n^ y^^/y^-1 

In the latter expression the second term in the r.h.s. is due to the breaking of the mass 
degeneracy. With the use of ( p.34| ) the energy ( |2.33| ) can be casted as 

E = ^ 'dy y +^ ''dy y . (2.35) 



Ji 27r72^9^ ^2/^-1 



However, notice that (|2.34| ) can only describe configurations for which the string is longer 
than 

Mq 

V? dy 



■'lim 



2T:mq J I y'^^/y'^ - 1 



(2.36) 



This value obtained for Uq = 2'Kmq. For L < Lum, the string does never reach values of U 
smaller than 27rmq, as can be seen in fig. 4. Therefore, the separation between the string 
ends and the energy are given simply by 

L^^r^ E^^r^,y^. (2.37) 



E 




Fig. 5: The energy in terms of the length for rriq = 1 and, from bottom to top, 
Mq = 1,2,3, l.h.s. For the Mq = 2 case we show the asymptotic (2.41) and (2.42). 
In the right panel we show the case of infinite heavy mass for i) the D3-D7 system 
(2.38), thin curve and ii) for the D4-D6 system (3.29) thick curve. 
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In fig. 5, we have plotted the energy as a function of the distance between the string 
ends. For large L, the effect of incorporating a heavier mass shifts the curves upwards but 
without modifying its characteristics. But for small L a novel effect w.r.t. the equal mass 
case appears: instead of a pure linear behaviour there is a slight positive curvature, a kind 
of parabola, that increases with the mass difference. Notice that in this case the slope, as 
approaching L —>■ , decreases indicating a smaller running of the /3- function. 

2.2.1. Mass Subtraction 

Notice that ( p.35| ) is not a renormalised expression, and therefore it will diverge as 
soon as the upper cutoffs are moved to the boundary at infinity. In that case, we shall 
deal with an infinitely heavy, degenerate pair of quarks. To render ( |2.35| ) finite in this 
case, we subtract the mass of the static quarks given by a couple of strings, one for 
each quark, stretching from the boundary to the horizon along the transverse coordinate 



Different shapes in the string providing the subtraction must correspond to different 
renormalisation schemes in field theory. We have checked, in the degenerate case, that 



the string configuration: we locate a pair of quarks, qq, on the D7-brane with a fixed 
distance among them, Lq. The dynamics is purely dictated by the motion of the qq system 
under the potential created by the stack of the D3-branes at the origin. Minimizing 
( pj.35| ) and ( pj.38D as a function of m^, E{Lq, niq), one can verify that the energy of the 
system is minimum when ( |2.35|) is used. In increasing niq to high values the difference 
between ( p.35|) and ( |2.38| ) is reduced until we reach the limit niq oo where the role is 
interchanged and the energy minimum is given by (|2.38| ). 

It is instructive to discuss some approximation to different scales: 

2.2.2. Long strings 

At low-energy we shall take the following set of assumptions 



Solving iteratively (always at leading order) one can find the turning point for a given 
string length 



U{a) = a, Xo = r [1. With this (12:351) reads 




(2.38) 




for finite masses gives the minimum energy for 



Length = L fixed and large ; niq , Mq ^ 2n . 




(2.39) 
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with C given in ( |2.21|) . For Mq — > we recover the symmetric case [15]. Corrections 
due to heavy masses are suppressed. 

Within the same approximation the energy gives 



E = Mq + mq + 



Uo 



Uo 
2n 



dy 




(2.40) 



Without performing the hmit Mq, rUq oo there is is no need of bare mass subtrac- 
tion as in [H,^ (see discussion above). Using the leading order in ( p.39|) we obtain 



E ^ Mq + ruq - 

1j 



1 - 



1 n^c^ 



12 (7rL)3 \m 



1 1 

+ 



^2.41) 



Notice that the the binding energy between quarks is attractive and corresponds at leading 
order to a system of two nonrelativistic particles. 



2.2.3. Short strings 

We shall consider now the case where the string is short. The limit L ^ is clearly 
achieved by taking Uq 0, ( p.37| ). Then, it is easy to prove that only keeping the first 
nontrivial contribution of the Taylor expansion of the energy in terms of the string ends 
separation, one has a quadratic behaviour 



E^M. 



67r2 



M3 



(2.42) 



contrariwise to any field theory expectation. 



3. QCD3+1: D4-D6 brane model 

A more realistic model, in the sense of being four-dimensional and not supersymmetric, 
can be obtained by placing a stack of Nc D4 branes at the origin and wrapping them on 
a 5"^ cycle. The latter are the source for fermion mass at one loop once antiperiodic b.c. 
are assumed []T^. As previously we can set an Nf L>6-probe branes displaced a distance 
in some subspace orthogonal to both brane sets. 

The field components at leading order in a' are 
i) The lOd background metric 

"^"'o = (I) ^ (""^^^ + "^"^^ + - + "^"^^ + ^ ^^^^ "^^^ 
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with f{U) = 1 — , = (irNcgs)^^^ a'^/^ and Uh denoting the horizon location. 
ii) A constant four-form field strength (in string units) 

^ R^LJ4, (3.2) 



with LJ4 been the volume form of the transverse S"^. 
Hi) Finally, the dilaton field 



e* = * I ^ ) . (3.3) 



3/4 



R J 

The relation with the gauge field theory parameters is given by 

a 



To be in the supergravity approximation the following requirements of the fields 
have to be fulfilled 

i) The smallness of the scalar field ( p.3|) . This restricts the reliability of the background 
to the infrared region. 

ii) The smallness of the curvature invariant and higher order derivatives. 

This two restrictions amounts to take the 't Hooft limit Qym-^c ^ 1- 

In addition, in order to describe QCD3_|_i one has to demand 
Hi) that the energy scale is such that its associated transverse coordinate U in the 
supergravity description corresponds to values where the background is indeed dual 
to QCD3_|_i rather than QCD4_|_i. 

The third restriction can be examined from the perspective of the quark-antiquark 
energy, which will be considered in section (3.3). In fact, it is only in the regime of 
large distances L between the quark and the antiquark that we reproduce the area law 
describing the feature of confinement of QCD. At small distances the 5"^ opens up and 
the background no longer describes QCD3_|_i but QCD4+1. If the position of the D6 
brane, where the quark-antiquark stay, is far away in the transverse coordinate, we shall 
deal with a Coulomb phase for the 4 + 1 field theory. Large distances means that the 
string of the Wilson loop goes deep into the bulk along the transverse coordinate until 
it bends in a minimum near the horizon. The bending must occur at some Uq such that 
^°~^>^ ^ 1. This is a third condition that must be met in order to stay in the region 
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where the background exhibits the area law behaviour for the Wilson loop. Note that this 
condition will guarantee that U ~ Uh Mkk, which is the regime where the theory is 



effectively four dimensional [|T9 



Finally, let us point out that the mass of the quark is given by the energy of a string 
stretching from the horizon to a flavour brane 



m, = — dU . (3.4) 

3.1. Spinning mesons 

We turn once more to the analysis of rotating strings. The initial point is (|2.3|) and 
the string conflguration is identical to ( |2.6D . The relevant part of the metric is 

+ (3.5, 

That leads to the action PO] 



^ ^ 2^ / W""' (c^^ 5^)- 

As in the D3-D7 case the variation of the action imposes that the strings end orthogonally 
to the D6-branes. 

In general the e.o.m derived from ( |3.6| ) are difficult, if not impossible, to be tackled 
in an analytic way, but as previously the solutions can simplify by choosing the proper 
gauge fixing conditions (k = 1). 

i) U{a) = a 

(3.7) 



a) r{a) = a 



R^uj^a 3 fU?-2U 



2 



3 /2f/3-f/3\ 



U^ + ———^ — 77^ = 0- (3-8) 



In addition we impose the b.c. = fixed , 9crL''{^g|o-=±f ^ oo , (z = 1, 2) 
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The invariance of (|3.6| ) with respect to t and 9 signals the energy and spin conser- 



vation 



E 



3/2 



,2 I R^U'' 



\ 1 - (a;r)2 



J 



2'Ka' / \R 



3/2 



„/2 I fl3t//2 

\| 1 - [ujrY 



(3.9) 



(3.10) 



For the rest of this section, we set 2'Ka' = 1 . 



3.1.1. Meson spectrum at large J 



Once again, the regime of large angular momentum corresponds to long strings 
rotating with small angular velocity u ^ Q . 



-1 -0.5 



Fig. 6: The profile for long strings, u) ^ Q . Equal mass case and different masses 
(l.h.s. and r.h.s. respectively). The dotted line denotes the horizon position, the 
horizontal full lines the flavour brane location and the dashed line the analytic 
string continuation. 



In fig. 6 we show a numerical solution for the string profile. This is approximately 
given by that of a static long string 

dur^iRUQ?^'^ ^ (3.11) 



which is characterised by almost straight strings. Notice that (|3.11D diverges in the limit 
Uq — > Uh- As previously Uq denotes the string turning point. 

The numerical solutions to (^.Tj), ( ^.8) ) are almost identical to those depicted in 
fig. 2 and fig. 3. Even possessing these numerical solution is not too much illuminating. 



As first noticed in [20|, an inspection of fig. 6 reveals that for slowly rotating string. 



J ^ yfg^Wc, a good approximation is to consider three parts of the string: two stretching 
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in the U direction from U^q and C/^g to Uh and the third horizontal piece where the string 
stretches in the field theory directions (r) at U = Uh- This allows to split approximately 
the integrals ( p.9|) , (|3.10| ) in three regions. The two vertical regions, correspond to r'^ ^ 

3/2 3/2 

Tj3_Tfi and the horizontal region {U^Uh) to r'2 > Therefore, the energy reads 

E = 



3 

= ( ^ ) (arcsin(wr^ ) + arcsin(wrMj) -\ , H , , (3.12) 



whereas the angular momentum can be expressed as 

J = 



3 

^"^ (^arcsin(wr^J - wr^^ y^l - wV^^ + arcsin(wrMj - wrM, y^l -wV^^J + 



+ , ■ (3.13) 

Notice that this is just a straightforward generalization of the results presented in ||2^ . 
Hereafter and r^, stand for the radius of rotation of the light and heavy quark 
respectively. 

Recall that not all of the parameters in the above expressions are independent. For 
given values of the masses, once u is fixed, the values of r^^, r^, follow from it as a 
consequence of the boundary conditions that must be imposed on the string. In order to 
compute this relation, we will follow the steps of |jl^ and define 

r{U)=r,t{U)+6r{U), (3.14) 

where rgt stands for the static configuration ( |3.11|) . Considering infinitesimal cu^, Sp we 
can expand ( |3.7| ) and get a first order differential equation in 5r'{U) 

{U'-U'o)iU'-Ul)6r''+^-^mu',-Ui)+^^^ =0, (3.15) 



* In order to see this, consider a = U, so both terms go to infinity but the first one dominates 
as the string profile at small u) must be similar to the static one, so ( |3.11| ) holds approximately. 
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which can be immediately integrated 

Sr'iU) ^ ,JL (-n^^^ r ^^^-m^^. 1 . (3.16) 



The integration constant has been encoded in f/ref? that we define as the value of U at 
which 5r'{U) = 0. The boundary conditions are that both ends of the string reach both 
fiavour branes orthogonally ( |2.12|) 

r'iU^ol) = r'sMl) + Sr'{ul;l) = . (3.17) 

As we are in the long string regime, Uq ~ Uh and the integral in (|3.16| ) is dominated by 
the vertical part of the string. Moreover, unless ^ oo, one can safely t/ref ~ f^o as the 
point where the derivative does not get modified must be in the horizontal region of the 
string. In the degenerate case = Mg, one has just f/ref = Uq because of the symmetry. 
The relevant integral is, then J^^*^ lia; ^ i-^^r^ (Uoe - Uh) ■ From (|3.17|) , one 

gets two conditionsf 

1 - ^^^li, = ^"^Tm^ rhq, 1 - uj^M^ = T-^ uj'^tm^ Mq , (3.18) 
where Tg is the tension of a string stretched aX U = Uh 



3 

Uh^^ 
R 



(3.19) 



and we have defined for each quark 



ft=^^^fF^- (3.20) 



Notice that for Ud6 ^ Uh, one has m ^ m (|3.4|), i.e. the mass of the quark, recovering 
the results of [^, which are the holographic realization of the hadron model of a rela- 



tivistic string rotating in flat space [21,^|. Identifying m, M with the mass of the quarks. 



( p.l8|) ensures the conservation of the center of mass motion for a rotating string with 



massive endpoints attached [21 



mqrmr.'^ Ts r r MqVMa^ T, 



'-^^^ - - - rl^y = ' : - -a/1 - rV^' ■ (3-21 



\ Notice, however, that the procedure to obtain the constraints (3.18) differs from that in 



In [20|, the string was cut into three pieces and (3.18) is the condition for a proper gluing 
of the three pieces. Here, the string was treated as a whole and (3.18) comes from the boundary 
condition of the ends of the string hitting the branes orthogonally. 
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This is certainly signaling that gluons inside the flux tube play a dynamical role 
and not only act as a mere geometric constraint, contributing to the total energy of the 



system. If not, momentum conservation would be just (consider r^^w^ ^ 1) 



(3.22) 



as in the supersymmetric case. In this nonrelativistic limit, the momentum conservation 
reads 

(3.23) 



"q' mq I 2 "^9 ""9 ' 2 

We expect a smooth transition between both behaviours, ( |3.22|) , and ( |3.23| ) 




Fig. 7: Chew-Prautschi plot obtained numerically for Uh = -2, Ujjg(i) = 2tt and 
Ujjq{2) = 2tv (squares), U^q(2) = 47r (triangles). The solid lines obtained using 
( 3.12| ), ( |3.13[ ), ( p.l8| ) are the extrapolation of these curves to the large-J region. 
For comparison in the r.h.s. we have plotted the experimental nucleon (full curve) 
1 23] and Pomeron (dashed) ||23] Regge trajectories. 



In fig. 7, we have depicted the Chew-Prautschi corresponding to this model. The 
curves giving the large-J limit are represented along with some points obtained numerically 
from (|377|), (|3.8|), (p^), ( |3.10| ) in the region of intermediate J. Unfortunately, we have 
not been able to perform a trustable numerical analysis of the very large-J region, because 
when the string has a nearly vertical-horizontal- vertical profile, numerical instabilities ap- 
pear. Despite the long extrapolation, the numerical results are robust and trustable in the 
intermediate range of J we considered. For comparison in the r.h.s., we have presented 
some experimental data. Notice that the qualitative similitude with the l.h.s. in the full 
region is remarkable. 



Given the above naive picture it is tempting to think in a very special limit of QCD, 
heavy quark effective theory, consisting of heavy hadrons composed by one heavy quark 
and one light one ||2^. The heavy quark momenta is given by Pq = itiqVq + , been 
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/c^ a residual momenta of order Aqcd ■, while the hadron total momenta is Pj^ = mHV^. 
In the heavy quark limit, the velocity superselection rule leads that the residual momenta 
is negligible and all the hadron momenta is carried by the heavy quark, Ph ~ Pq Il26[] . 



This setup is certainly susceptible to be obtained by the above model using (|3.23|) via the 
identification 

- m,T^^ + ^rl^ ^ 0{Aqcd) • (3.24) 



Notice that ( |3.24| ) is suggesting that for a given mass the average position of the light-quark 



with respect to the heavy one is not arbitrary. 



3.1.2. Intermediate J strings 



In the region of intermediate J, the holographic model differs from the string ro- 
tating in fiat space [p0| , |2l| , p2l because the string profile cannot be approximated any more 



by vertical and horizontal regions. It may be regarded as a string of nonconstant tension 



rotating in fiat space. Similar models where considered in in the search of more 

realistic models. Here, the deviations from the constant tension are given uniquely by the 
dynamics of the rotating string of the ten-dimensional theory. In fig. 7, the points obtained 
numerically show the deviation of the formulae computed in the large angular momentum 
approximation when J goes to intermediate values. 

It is tempting to consider these deviations as the holographic realization of the 
corresponding ones in QCD. However, one should keep in mind that out of the large- J 
limit, quantum corrections to the semi-classical model must also be taken into account. 
These kind of corrections where computed in for a static configuration. Our setup is 



much more involved so it seems hopeless to perform the analogous computation. 

3.1.3. Small J strings 

For completeness, let us consider the J — limit, i.e. rotating short strings. 
However, when the string is short, it probes the KK direction, so one cannot think of this 
setup as the holographic small- J dual of QCD3_|_i. By taking r ^ in ( ^77|) one finds 



r —r 5 ^r = 0. 3.25 

We have not been able to find an explicit solution for this equation. However, a numerical 
computation shows that the qualitative features are similar to those of section 2.1.1. for 
the supersymmetric case (see fig. 2 ). 

3.2. Quark- antiquark energy 
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Fig. 8: The quark-antiquark potential as a function of the distance for flavoured 
branes with, from bottom to top, ruq = 1 and Mq = 1,2,3. We also show the 
asymptotic limits (3.30), (3.32) and (3.34). The r.h.s. shows for comparison the 
unquenched lattice results p9[ . The squares denote the ground state static poten- 
tial while triangles stand for the first excited state. 



In this section we change tack, retreating from the semi-classical limit of spinning 
strings to make qualitative statements about the behavior of the static quark potential. 

The D6 brane configuration has its world- volume along xo,xi,X2, xs and along two 
coordinates of the S^. The setup is realized by the string configuration 

U = U{a), Xo = r, Xs = a , (3.26) 

where we have stuck to the static-gauge. The expressions for the energy and the quark-anti 
quark separation are basically the same as those obtained in but with an ultraviolet 
cutoff, equal to the positions of the quarks on the D6-branes 

yd) y(2) 

i?3/2 I ^3/2 I 

c/o/ A Viy^ - 1) (y^ - ul'^J-^ Viy"" - 1) (y^ - v^) 

yd) y(2) 

J, Up f^, , ^0 f^, 

h = — / ay — = H / (n ay — = , ri = — . 

TT A V(2/3 _ 1) (^3 _ ^3) 27r y-ki ^(^3 _ 1) (^3 _ ^3) 

(3.28) 

In fig. 8 we display the relation energy versus distance in the qq system for a given pair 
^dI ' '^^^ gross features for the ground state, convexity, match those of QCD 



but with a relevant salient point: the existence of a linear behavior at large distances (see 
discussion below). As in the D3-D7 case the mass splitting only shifts upwards the curve 
E vs. L in the large L region without introducing any other change. At high-energy the 
nondegenerate case shows once more a parabolic behaviour. In fact the profile of the l.h.s. 
plot is identical to the one in the supersymmetric case. 
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3.2.1. Mass Subtraction 

As in the supersymmetric case the energy ( p.28|) is not a renormahsed expression 
once the integral upper cutoff go to the boundary, and needs subtraction. With the same 
procedure as previously ( |3.28| ) reads 

- '"^ *Wr^. (3.29) 




Which, incidentally, does not match with the results of [3C]. 



As it is shown in [jT^ , and unlike in the conformal case, now the D6-brane is not 
located at a fixed point in the transverse U variable, but rather bends outwards in the 
region where it is closest to the D4 stack of branes. We have checked that the qq system 
minimizes its total energy when located at the place in the D6-brane that minimizes its 



distance with respect to the D4. This result coincides with the conclusions in [jT^. In 
addition, as previously, in order to obtain a light, finite qq system we have to make use of 
instead of (Q). 

3.2.2. Long and short strings 

We turn to discuss the behaviour of long and short strings. Notice, however, that 
although the relative distance between the location of the qq pair on the D6-brane and the 
horizon can be small, this by no means leads to the conclusion that the distance between 
q and q must be short as one can indeed verify from ( |3.27| ) by approaching the horizon at 
fixed C/jJg, i.e. y — > 1. In fact for ry, y ^ 1 ( |3.27| ) develops a single pole that can not be 
cured, and hence the string approaches an infinite length. 

Long strings 

A feature matching the expected QCD behavior is found at low-energy. There the qq 
system has been hadronised and from the macroscopic point of view we deal with hadrons. 
The slope of these strings is constant. This is consistent with the Regge behaviour shown 
for rotating long strings in section 3.1.1. In order to obtain long strings at a given finite 
distance for U^^^ ^ (z = 1,2) the string turning point should approach the horizon. In 
that case ry = 1 + e and the relation between the energy and the length along the brane 
takes the form 

^-i(^)^'^^-2^Ki + f/S-2t/.). (3.30) 

Notice that the inclusion of heavy masses shifts the intercept but does not modify the 
slope, thus we still deal a Regge behaviour as in the case of static quarks [0. The string 



tension is the same as the one obtained in [30|, and matches the expression obtained from 
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the energy formula for rotating strings . The reason is that at large distances the main 
contribution to the integrals ( |3.27| ), ( |3.28| ) comes from the vicinity of the horizon, so the 
result is insensitive to the presence of the cutoff at Uj-)^. 

Short strings 

As in the supersymmetric case, when the QCD string is very short, its holographic 
dual is given by a string directly stretching between both flavour branes, without bending 
(see fig. 4). Then, the expressions for the distance between the ends of the string and its 
energy are just 
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Clearly, the L corresponds to taking the integration limit Uq ^ (notice that this 
makes sense since in this case Uq does not have the physical meaning of the string turning 
point). Then, as in the supersymmetric case, one can check that, at leading order, the 
energy has a parabolic dependence on the distance 

E = Ci+ C2L2 + . . . , (3.32) 

where the constants are given by 

1 C6 X2 1 

(3.33) 

For the sake of completeness, let us mention that, in the degenerate case, this reasoning 
is clearly not valid and the string must stretch with bending. Then, the leading order 
expansion leads to a constant slope behaviour for the energy 

E = — (^Y\ + .... (3.34) 




27r V ^ . 

This shape has been quite long discussed in QCD, see for instance [^. With the range 
of energy explored in the actual lattice data this effect can be interpreted as a renormalon 
artefact and has no further physical consequences. In the present case we can understand 
it in an energy basis: in the coordinates (|3.1|) E U . As we increase the U coordinate 
also the energy is increased. At certain energy, E ~ Mkk ~ §'^172" there is an opening of 

the dimension on the , thus we deal with QCD4+1 instead of QCD3+1. 

To substantiate this claim in an analytic way, we have performed the limit t/^g 
00, and renormalised the energy with the subtraction (|3.29| ). We can observe that for small 
distances between the quarks, the energy behaves as 

E^-7]U^/^^+ct., 
thus showing a Coulomb phase for QC 04^1. 
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3.3. Quark- antiquark potential screening 

It is instructive to gain some insight on the effects of including the dynamical quarks, 
mimicked for the D6-brane. In fig. 5, r.h.s., we have plotted (thick curve) the quark- 
antiquark energy obtained by taking Ud& ^ oo and subtracting the bare mass ( p.29| ) . As 
one can see, already keeping the efi'ective field theory in (3+l)d, the effect of introducing 
dynamical quarks is quite noticeable and with similar pattern to the supersymmetric case. 
This change of behaviour can be understood as a screening effect: while in a pure glue 
dynamic system composed of two infinite external heavy-quarks the interaction can be only 
mediated via gluons, in a system with fundamental matter the gluon propagator itself is 
dressed with quarks. This fact weakness the interaction between the valence- quarks and 
the potential wall gets broadest. To manifest this effect we have fitted to fig. 8, inside the 
(3-M)d region, two phenomenological models for the case rUg = Mq = 1. The first model 
is given by the screened potential pc 



V{r) = (-- + ar) , (3.35) 

V r / fir 

and the second by the well-known Cornell potential 

V{r) = Vo + br--. (3.36) 
r 

In this case the model realised with D6-branes behaves most as a screened system at short 
distances. Although the relative different between potentials is so small that no conclusive 
statement on the behaviour of the system can be obtained. Similar conclusions are achieved 



by lattice studies i34|,|35 



As a concluding remark on this section we mention that for both, the singlet and 
the first excited state, see the UV region in fig. 8, the beta function for the D4-D6 system 
has an UV fixed point. This is in the singlet case is in clear contrast with our actual 
knowledge on perturbative QCD. 

3.4. 't Hooft line: monopole-monopole interaction 

A complementary way of looking at confinement is via the magnetic dual of the 
quark, the monopole. The configuration realizing a magnetic monopole from (|3.1| ) is that 
of a D2-brane ending on the initial D4-brane. The monopole is obtained by wrapping 
the D2-brane on the compactified |3^,^. With our conventions, the mass of a free 



monopole is 

mr,=KiU^06)-Uh). (3.37) 

where we have defined the constant K = AitTjj2R^^^Uj^ ^^"^ 9^^ • The potential between the 
monopoles is obtained by placing two monopoles at distance L in the xi-axis. Consider a 
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D2-brane parametrised in term of (t, U{xi), 62) ■ The induced metric on the brane is then 
3| 



9u = - {^"^ , 9UU = {^^Y^ f{U)-\l + f{U)x\{Uf 



Inserting it in the D2 action 



m— m 



Td2 j dxid02e ^/gI^guuge^ , (3.38) 



we can extract the distance between the ends of the flux-tube and the energy stored in it 
(for simphcity, we consider monopoles with degenerate mass) 

ul'^ Ji Viy'-^)(y'-v') ' 




E^KUol dy\l'-^^, (3.40) 
where 77 is defined as in (|3.28D. 



For small distance there is a configuration of the flux tube between the monopoles 
for which the energy (|3.40|) is less than that of two free monopoles (|3.37|) . As we increase 



the distance between the monopoles the energy stored in the tube grows up to certain Lcrit 
where it reaches 2mm signaling the breaking of the flux tube. From these the most stable 
configuration corresponds to two flux tubes stretching from each of the monopoles in the 
D6-brane up to the horizon. This leads to a constant energy density in accordance with 
field theory expectations . 



Notice that for probe branes near the horizon, light-quark masses, the breaking 
takes place at lower values of Lcrit, while for probe branes supporting heavy quarks the 
breaking occurs at higher ones. 



4. Summary 

By considering open rotating strings attached to probe branes we have evaluated 
semi- classically some properties of a system mimicking excited mesons formed by different 
constituents masses. We focused on two basic aspects: i) the Chew-Frautschi plot and 
u) the qq system energy. The first quantity is obtained by introducing open rotating 
strings with ends attached to an embedded probe brane an extending into a transverse 
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coordinate. While the second is a static quantity. In the actual range of energies the QCD 
static potential has the form 

= -- +cr. 

Adding light flavours we move away of this limit and new effects, as the spin, play an 
essential role. These potentials are reliable to study heavy quark interaction, but when 
light quarks are involved, relativistic effects are important and even the concept of inter- 
quark potential remains questioned. 

Bearing this lack in mind we construct a system composed by two different, tuneable 
masses by adding two probe branes and separating them in a plane orthogonal to the initial 
brane stack world- volume. 

The first case of backgrounds we treated is a supersymmetric setup, the D3-D7 



brane system, generalising the extremal case |15]. At low-energy neither the J — E relation 



nor the E — L show a qualitative difference with respect to the degenerate case but at 
moderate-high energy differences are substantial. Notice that at qualitative level, and 
without any apparent reason, the J — E and E — L plots show the same shape, see fig. 1, 
fig. 5. 

The second studied background is non-supersymmetric and is obtained considering 
a thermal cycle in lOd. The gross features of the Chew-Frautschi plot for hadronised sys- 
tems resembles that of heavy-quarkonia in QCD. The E — L relation presents a Regge type 
behaviour for large separations, matching the field theory expectations. At high-energy 
there is however also a Regge behaviour in the case of degenerate masses. This result can 
be interpreted, granus salis, as the possible presence of a screening effect as one can see 
by fitting the qq energy to two different potentials. There is however a rather surprising 
quantitative agreement with lattice findings if the singlet and 1st excited states are con- 
sidered. The basic tenets of the model can be tested by lattice Monte-Carlo simulations of 
higher excited states. The last remark we have made on confinement concerns the picture 
of a monopole-monopole condensation, which fulfills the field theory expectations. 

We stress that the all the mentioned features are independent of having imple- 
mented supersymmetry in the theory. In this respect we believe that the aforementioned 
characteristics of the Chew-Frautschi plot or in the quark-antiquark potential in the D4-D8 
||40|| model must led to the same qualitative results as those in the D4-D6. Nevertheless it 
would be interesting to see if the existence, at large-energy, of a Regge behaviour persists 
in this model. 
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